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Abstract. Non-minimally coupled scalar field models suffer of unstable growing modes at the linear perturbation level. The 
nature of these instabilities depends on the dynamical state of the scalar field. In particular in systems which admit adiabatic 
solutions, large scale instabilities are suppressed by the slow-roll dynamics of the field. Here we review these results and 
present a preliminary likelihood data analysis suggesting that along adiabatic solutions coupled models with coupling of order 
of gravitational strength can provide viable cosmological scenarios satisfying constraints from SN la, CMB and large scale 
structure data. 
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INTRODUCTION 

The possibility of a direct coupling between a quintessence-like scalar field and the various matter components has 
been extensively studied in a vast literature (see e.g. [1]). Non-minimally coupled scalars appear in various theoretical 
scenarios which attempt to describe fundamental interactions at energies beyond that of the Standard Model of particle 
physics (see e.g. [2]). Their application to cosmology and the unsolved problem of dark energy in the universe has 
suggested a number of interesting features, most importantly the solution of the so called coincidence problem. The 
presence of the scalar interaction is not incompatible with existing constraints on the violation of the Equivalence 
Principle. Ingenuous mechanisms (which may differ from one model to another) guarantee that standard General 
Relativity is recovered on Solar System scales (e.g. [3, 4]), and leave distinctive signatures on the structure formation 
process, eventually contributing to some of the still not understood phenomena in the context of Cold Dark Matter 
paradigm [5]. Consequently the distribution of structures, at least on those scales where observations have provided 
accurate measurements, is a key test that such models have to pass. Nevertheless a number of works have indicated 
that coupled scalar field models may suffer of large scale instabilities at the linear perturbation level. This was initially 
pointed out in some specific realizations [6] and recently discussed in more general setups [7, 8]. The claim has been 
particularly emphasized on scenarios characterized by the existence of the so called "adiabatic" regime, such as the 
Chameleon model [7]. In the light of these results coupled models seem to be unrealistic cosmological scenarios. 
However as we have shown in [9] the instabilities are not generic, rather they are strongly dependent on the dynamical 
state of the scalar field. More importantly along "adiabatic" solutions and for natural values of the coupling constant 
such instabilities are innocuous. Here we will briefly summarize the main results of [9] to which we refer the reader 
for a more detailed discussion. We will also present the results of a preliminary likelihood data analysis to test the 
viability of these models against current cosmological observations. 



PERTURBATIONS IN COUPLED SCALAR FIELD-DARK MATTER MODELS 

Let us consider a scalar field 0, with potential y(0), coupled to matter particles through a Yukawa couphng of the 
form f{(p /Mpi)\i/\fr, where / is the coupling function and y/ is the Dirac spinor associated with the matter particle 
(Mpi — 1 / vSttG with G the Newton constant). For simplicity let us consider the case in which the scalar field is 
coupled to dark matter only, thus Equivalence Principle constraints are immediately satisfied. This is not a restrictive 
assumption since our results can be extended also to models with couplings to all matter components provided the 
existence of an "adiabatic" regime. Because of the coupling the energy-momentum tensor of each component of 
the system is not conserved. It is only the total energy-momentum tensor that satisfies the conservation equation: 
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0. Now let us consider a coupling function of dilatonic type, f{<p) = exp{^<p/Mpi), with j3 
the dimensionless coupling constant, from the above conservation condition in a flat Friedmann-Lemaitre-Robertson- 
Walker background (ds^ = —dt^ + a{t)^dx^) we obtain the evolution equations: 



Pdm + ^HpoM = j3 pDM , 

<^ + 3H^ + V^ = -PpDM, 
with the Hubble rate given by = {fj^ = ^ [poM + <^^/2 + V(<|>)] . The solution to Eq. (1) reads as 
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where 0o is the present scalar field value. From Eq. (1) and Eq. (2) we may notice that for positive values of the 
coupling constant j3, the interaction transfers energy from the 0-field to the dark matter particles, with the scalar field 
evolving in an effective potential 
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characterized by a minimum 
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Given the above background equations, the evolution of linear density fluctuations can be studied by perturbing 
the Einstein equations and the conservation of the total energy momentum tensor about a linearly perturbed FLRW 
background. However in order to gain some intuitive insight on the behavior of the perturbations on the large scales 
and perform a simple stability analysis, it is convenient to consider the interacting scalar field-dark matter system 
as an effective single fluid with energy density pj = ^'^/Z + V ((^) + Pdm, pressure pr = ^^/2 — and whose 
perturbations are uniquely characterized by an adiabatic sound speed, CaT = \/ Pt/Pt, and the rest frame sound speed, 
CsT = \/dpT/5p,. 

In synchronous gauge the perturbation equations reads as 
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where 5t = dpj / Pt and Qj is the shear velocity perturbation of the fluid. For a barotropic component with a constant 
equation of state (e.g. matter, radiation) c] = c^ = w. This is not the case for a generic fluid (e.g. scalar field), for this 
reason we may expect the effective unified fluid to be non-barotropic, (i.e. c^j. ^ c\j ^ wt). In terms of the scalar fleld 
and dark matter perturbation variables we have 



3Hf- + (l>[2V,^,+l3poM] 
3H^^ + 3HpDM ' 

(j)5^ + V^^5^+PDM Sdm 
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These relations provide us with a simple way of determining the stability of the perturbations in the coupled system, 
for example in a given background regime instabilities may develop if these sound speeds acquire sufflciently negative 
values. 
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FIGURE 1. Scalar field effective potential at z = 0, 3, 10 and 10^, the dashed line shows the position of the minimum as function 
of the redshift. 



An attractor solution of the background homogeneous system consists of the field seating at the minimum of the 
potential, and drifting in time according to Eq. (5). This is usually referred as "adiabatic" regime. It has been shown 
in [10] that along this solution the field slow-rolls, thus it has a negligible kinetic energy. In particular for a power law 
potential, V{<j)) 'x<j)~", the evolution of the scalar field given by the condition Eq. (5) reads as: 



which depends on both the slope a and the coupling j8. Equation (10) is a non-linear algebraic equation which can 
be solved numerically through standard bisection methods. The presence of the minimum distinguishes two different 
sets of initial conditions: < ^^l^ (small field) or (j)i„i > (j)^}^ (large field). For small field values, (j) evolves over the 
inverse power-law part of the effective potential, where it minimizes the potential by slow-rolling as shown in [10]. In 
contrast for initially large field values, (j) rolls towards the minimum along the steep exponential part of the effective 
potential. Thus it rapidly acquires kinetic energy which subsequently dissipates through large high-frequency damped 
oscillations around the minimum. The growth of the linear perturbations in the coupled scalar field-dark matter system 
is significantly different in these two regimes. 



Let us evaluate the adiabatic and rest frame sound speeds along the adiabatic solution respectively. Substituting 
Eq. (5) in Eq. (8) and neglecting the term proportional to the kinetic energy of the scalar field (due to the slow-roll 
condition) we have 
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Adiabatic Regime: slow-roll suppresion of instabilities 
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FIGURE 2. Upper left panel: evolution of the scalar field equation of state and effective unified fluid equation of state wj'. 
Right upper panel: evolution of the scalar field velocity with respect to the Hubble rate; Lower left panel: redshift evolution of the 
adiabatic sound speed c^j and propagation of pressure perturbations c^^ ; Right lower panel: Linear growth factor of the dark matter 
density contrast a.tk= 10~^, 10"'^ and 0.1 Mpc~^ . 



since ^ > it then follows that c^j < 0, implying that adiabatic instabilities may indeed develop. However we should 
remark that during the adiabatic regime the field is slow-rolling (i.e. 3H(j) « 0), hence the term (j) /3H can be negligibly 
small compared to j3, such that c^j- w 0^, thus leading to a stable growth of the large scale perturbations. In fact let us 
suppose that c^j. = —10^^, the instability will affect modes k > 10^, but these correspond to very small scales which 
are already in the non-linear regime and for which the linear perturbation theory does not apply any longer. In contrast 
large scale instabilities will occur if the coupling assumes extremely large values, j3 ^ 3H/^. This is consistent 
with the conclusions of [7], where the authors have shown that during the adiabatic regime perturbations suffer of 
instabilities provided that J3 » 1. However such situation would be extremely unnatural introducing a large hierachy 
problem in the gravitational sector since it would implying having a scalar fifth-force which is (1 + 2j3^) greater than 
gravitational strength. Guided by naturalness considerations one might expect that the dimensionless coupUng constant 
is of order of unity. Let us now evaluate the sound speed in the total effective fluid rest frame, Eq. (9) we have 

c?r = (12) 

^ P 5<l> 

assuming that the scalar field is nearly homogeneous, 5^ <C 5dm (in Planck units), we have c^j « jid^/SoM, and for 
j3 PS ^(1) this implies c^.j « 0. In other words if the scalar field fluctuations are small with respect to the dark matter 
density contrast, then the coupled system behaves has a single adiabatic inhomogeneous fluid (c^j. « c^j. w 0). These 
results are supported by the numerical study of the perturbation equations for the individual components of the system 
as suimnarized in Fig. 2. 



Non-Adiabadlsc Regime: large field oscillations and onset of instabilities 



For initially large field values, <j) rolls along the steep exponential part of the effective potential. Its evolution is 
therefore dominated by the kinetic energy and the field behaves as a stiff fluid (w^ = 1, as can be noticed in the left 




FIGURE 3. Upper left panel: evolution of the scalar field equation of state w,j, ; Right upper panel: evolution of the scalar field; 
Lower left panel: evolution of the field fluctuations 50^ at = 10~^ , 10~^ and 0. 1 Mpc~^ respectively; Right lower panel: evolution 
of dark matter density for A:- values as in the case of d(j>k- 

upper panel of Fig. 3). Then as the field reaches the minimum, the kinetic energy is damped away through a series 
of high-frequency oscillations. During this oscillatory regime, which is similar to that of the inflaton in the reheating 
phase, the scalar field perturbations are unstable and exponentially amplified by the background-field oscillations 
provided that their frequency increases as their amplitude diminishes [11]. This is indeed the case as shown in Fig. 3, 
where we plot the evolution of (right upper panel), 50^^ and Sdm for three different wave-numbers, k = 10^^, 10^^ 
andlO-^We can see there that an instabihty occurs roughly at the same time of the first minimum-crossing oscillation, 
then followed by a second stage of exponential growth at the beginning of the second oscillation. Such unstable modes 
are similar to those found in [8], in fact by averaging over periods of time larger than the characteristic time of the 
oscillations, the scalar field behaves effectively as a dark energy fluid with a constant equation of state w, as the case 
considered in [8]. 

CONSTRAINTS FROM SN lA, CMB AND LSS: A PRELIMINARY ANALYSIS 

Coupled models have been tested against cosmological observations in various works [12, 13, 14]. The main con- 
clusion of these analysis is that current measurements of the CMB anisotropy power spectra and the matter power 
spectrum from galaxy surveys constrain the coupling constant to be j3 < 0.01 — 0.1 depending on the specific model 
reaUzation. However none of these works have considered the case in which the scalar field evolves in the adiabatic 
regime. To this purpose we used the lastest SN la-UNION dataset compilation [15], WMAP-5 years data [16] and 
matter power spectrum measurements from SDSS-5 data release [17] to test the viabihty of a non-minimally coupled 
scalar field model with power law potential in the adiabatic regime. For this purpose we have specifically set the field 
evolution to satisfy Eq. (10), implemented the perturbation equations in a properly modified version of the CMBFAST 
code [18] and run a Markov Chain Monte Carlo likelihood evaluation. For simplicity we assume a flat universe and 
fix J3 = 1 and a = 0.2, while let all other parameters to vary. Here we simply aimed to test whether an adiabatic 
solution can provide compatible fit to the data, and we leave to a future study a more detailed analysis of the full model 
parameter space, including a and j3. The marginalized ID likelihood are shown in Fig. 4 and the best fit value and la 
errors are: CIdm = 0.222 ± 0.025, Q.bh'^ = 0.02226 ± 0.00068, = 0.75 ± 0.03, T = 0.072 ± 0.017, = 0.89 ± 0.02, 
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FIGURE 4. MarginaUzed ID likelihoods. 



As = 0.75 ± 0.03 and b = 1 .44 ± 0.49. These constraints are consistent with those derived for LCDM cosmologies, 
the total is close to that of the vanilla LCDM model such the two scenarios are statistically indistinguishable. As 
shown in [10] differences due to the scalar interaction may arise on the small scale clustering of dark matter. Over- 
all this preliminary analysis suggests that adiabatic coupled models with natural coupling are consistent with current 
cosmological observations. 
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